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Abstract 


A  three-dimensional  integral  boundary-layer  approach  is 
developed  and  coupled  with  the  streamline  method  for 
theoretically  determining  the  vortex-type  flow  separation.  The 
governing  equations  are  solved  in  a  streamwise  system  using  a 
power-law  profile  for  the  streamwise  flow  and  the  Mager  pro¬ 
file  for  the  crossflow  .  The  reduced  ordinary  system  is  then 
coupled  with  the  streamline  equation  and  integrated  using  the 
fourth-order  Runge-Kutta  scheme.  Crossflow  derivatives  are 
evaluated  and  accounted  for  during  the  integration.  A  prolate 
spheroid  at  incidence  in  an  incompressible  turbulent  flow  is 
considered  as  a  test  case.  Good  comparison  between  the  theory 
and  (he  experiment  has  been  observed  for  the  case  of 
a  =  10  degrees.  For  the  case  of  high  incidence,  the  method 
predicts  a  qualitative  trend  but  deviates  quantitatively  due  to 
the  large  crossflow  involved. 


Nomenclature 


0,s 


crossflow  momentum  thickness 


=  crossflow  momentum  thickness 


=  crossflow  momentum  thickness 


9  =  streamline  angle 

A  =  tan  fi 

p  =  viscosity 

i.  n.  i  =  streamwise  coordinates 
p  =  density 

r  =  shearing  stress 
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h, 


h. 


*2 

M 

P 

Re 
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Li 

u,  v 

V. 
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x,  y,  z 
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d| 


major  and  minor  axes  of  an  ellipsoid 

entrainment  coefficient,  defined  in  Eq.  (9) 

skin  friction  coefficient 

local  body  radial  distance  from  the  centerline 

shape  factor,  d,/0u 

(H  -  l)/2 

metric  coefficients  for  coordinates  i,  n 

curvature  parameter,  _L_ 

h,hs 

I  3h, 

curvature  parameter,  — ~ 

h,h,  dt, 

Mach  number 

static  pressure 

Reynolds  number,  pV^  a/p 

distance  along  and  normal  to  a  streamline 

velocity  at  edge  of  boundary  layer 

velocity  components  in  streamwise  coordinates 

freestream  velocity 

body-oriented  orthogonal  coordinates 
Cartesian  coordinates 
angle  of  attack 

wall  shear  angle  between  external  inviscid  streamline 

and  corresponding  limiting  streamline 

ratio  of  specific  heats 

boundary  layer  thickness 

streamwise  displacement  thickness 

r  m)« 

crosswise  displacement  thickness 


streamwise  momentum  thickness 


The  research  was  supported  by  the  Naval  Air  Systems 
Command  (AIR-3 10D,  AIRTASK  WR023-02)  and  by  the 
Independent  Research  Program  at  the  David  Taylor  Naval 
Ship  RAD  Center. 

•Research  Aerospace  Engineer,  Aviation  and  Surface  Effects 
Department.  Associate  Fellow  AIAA. 


Subscripts 
w  =  wall 

°°  =  freestream 

l.  n  =  A.  r)  direction 


Introduction 


In  three-dimensional  flows,  the  envelope  of  limiting 
streamlines  has  been  used  as  the  separation  line.  Thus  defined, 
the  determination  of  the  separation  line  has  been  found  to  be 
very  difficult  analytically  even  with  the  most  sophisticated 
numerical  techniques  available.  For  most  flows  of  practical  in¬ 
terest,  however,  the  separation  is  of  the  vortex  type  where  vor¬ 
tices  are  created  as  a  result  of  the  convergence  of  viscous 
streamlines  above  the  line  of  separation.  The  line  of  separation 
can  be  properly  determined  by  the  envelope  of  viscous 
streamlines  in  the  boundary  layer.1 

The  problem,  therefore,  is  to  determine  the  realistic- 
streamline  pattern  based  on  realistic  pressure  or  shear  stress, 
either  by  experimental  measurements  or  by  means  of  viscous- 
inviscid  interactions.  An  attempt  at  using  pure  inviscid  pressure 
along  with  a  simple  friction  model  has  had  limited  success.-  A 
more  general  theoretical  method  is  yet  to  be  developed.  The 
purpose  of  this  paper  is  to  describe  a  theoretical  procedure  for 
determining  the  vortex-type  flow  separation  using  a  three- 
dimensional  integral  boundary-layer  solution  coupled  with  the 
streamline  method.1  As  opposed  to  the  usual  boundary-layer 
approach,  this  method  deals  with  the  convergence  of 
streamlines  raiher  than  with  the  vanishing  of  skin  friction.  A 
prolate  spheroid  at  incidence  in  an  incompressible  turbulent 
flow  is  considered  as  a  test  case. 

Calculation  of  Viscous  Streamlines 

To  simulate  the  physical  flow  so  that  the  flow  separation 
can  be  detected,  it  is  necessary  to  consider  the  streamlines  in¬ 
side  the  boundary  layer:  sec  Fig  I.  It  is  assumed  that  the 
envelope  of  converging  viscous  streamlines  is  above,  but  close 
to,  that  of  the  converging  wall-limiting  streamlines  Therefore, 
the  line  of  separation  can  be  determined  by  the  envelope  of 
converging  steamline  equations  derived  earlier:1-  - 
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where  conventional  symbols  have  been  used  with  geometric 
variables  illustrated  in  Fig.  2.  For  incompressible  flows,  the 
potential  flow  solution1  can  be  used  tor  the  pressure  gradients 
3P  3x  and  3P  3$  in  Eq.  (3).  The  problem,  therefore,  is  to 
provide  proper  values  for  the  shear  stress  derivative  3r  34. 

Integral  Solution  to  Three-Dimensional 
Turbulent  Boundary  Layers 

A  simple  analytic  model  for  the  shear  stress  derivative  has 
shown  some  promise.-  However,  in  more  general  cases,  where 
the  flow  separation  is  predominently  of  the  vortex  type  a 
more  rigorous  theoretical  form  is  in  order.  While  it  i.  not 
straightforward  to  derive  a  universal  theoretical  friction  model, 
it  is  possible,  in  most  eases,  to  provide  the  required  friction 
values  numerically.  In  so  doing,  the  use  of  the  integral  method 
seems  to  be  appropriate  because  of  its  compactness  and  com¬ 
patibility  with  the  streamline  approach.  The  (longitudinal)  in¬ 
tegration  can  be  performed  along  a  streamline. 
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To  close  the  system,  streamvvise  and  crossflow  velocity  pro¬ 
files,  auxilliary  relations  for  the  skin  friction,  and  flow  entrain¬ 
ment  are  needed.  The  streamvvise  velocity  profile  is  assumed  to 
be  the  power-law  form  for  two-dimensional  flow;  that  is, 


HD" 


where  z  =  4-d,  h  =  (H  -  l)/2,  and  H  is  the  shape  factor.  The 
crossflow  velocity  profile  takes  the  form  suggested  by  Mager.4 
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Fig.  1  -  Streamlines  in  a  boundary  layer 
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Fig.  2  -  C  oordinate  systems 


where  A  =  tan  ft.  The  boundary-layer  solution  procedure  in  the 
present  work  is  similar  to  that  of  Smith/  but  differs  in  several 
aspects  which  will  appear  later. 

The  entrainment  coefficient  C.  is  given  as  a  function  of 
the  shape  factor  H  by  Green,6  witn  the  constant  term  being 
modified  from  0.022  to  0.032: 

C  ,  =  0.025H  -  0.032  (9) 

The  skin  friction  coefficient  is  related  to  the  shape  factor  H 
and  the  Reynolds  number  based  on  the  streamvvise  momentum 
thickness  Re,,  .6 
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where  Cfo  and  H0  are  functions  of  Re0| 
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and  the  Re0  and  C,  are  defined: 
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The  expressions  for  the  curvature  parameters  k,  and  K, 
are  derived  in  terms  of  quantities  along  a  streamline.  For  lv,, 
the  geometric  relationship  between  two  surface  coordinate 
systems  (x,  f)  and  (4.  r\)  is  used: 
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Taking  cross  differentiation  and  letting 


h2  cos  0 
f 


result  in 
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where  the  sum  of  the  three  terms  in  parentheses  is  zero,7  and 
the  inviscid  form  of  Eq.  (3)  has  been  used. 

To  find  K2,  the  continuity  equation  is  written  for  the  exter¬ 
nal  flow  in  both  the  x,  f  coordinate  system  and  the  4.  n  system. 
By  equating,  there  results 
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Numerical  Integration 


Equations  (4)  through  (6),  with  the  aid  of  Eqs.  (7)  through 
(12),  allow  solutions  for  three  basic  unknowns:  the  shape 
factor  H,  the  streamwise  displacement  ©,,,  and  the  tangent  of 
the  angle  between  the  external  inviscid  streamline  and  the  cor¬ 
responding  limiting  streamline  A.  These  equations  are  recast  in¬ 
to  the  form  of  streamwise  derivatives: 


DH 

Ds 

=  G, 

(13) 

D0l! 

Ds 

= 

(14) 

where  1)  Ds  -  (1  h,)0  J4,  and  Cl,.  Ci-.  and  G.  contain: 

(a)  geometric  properties  and  outer  edge  How  properties 

(b)  variables  H,  0n,  A,  etc.  and 

(c)  crossflow  derivatives  of  H,  0, , ,  and  A 

The  detailed  expressions  for  G,,  G2,  and  G,  are  presented  in 
the  Appendix. 

With  the  viscous  solution  in  hand,  the  crossflow  shear 
derivative  required  by  Eq.  (3)  is  obtained  by  differentiating 
Eq.  (8)  with  the  aid  of  the  following  relationships: 
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To  eliminate  the  unknown  eddy  viscosity  t  and  its 
derivative  3t/3z,  the  t  value  is  correlated  wiih  the  streamwise 
shear  stress  such  that 
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It  is  further  assumed  that  the  streamwise  shear  stress  varies 
linearly  across  the  boundary  layer: 

T{  =  V<*  _z)  =  y(?u2cf2  (1  -z)  1 

which  yields 

=  _  -2a  , 

yM2P  34  2d 

With  the  aid  of  relations  (18)  through  (21),  along  with 
basic  profiles  (7)  and  (8),  one  obtains 
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where  /  is  set  to  d|  d  for  actual  calculations.  The  choice, 
although  arbitrary,  is  considered  to  be  the  most  appropriate  in 
capturing  the  overall  property  of  the  boundary  layer. 

Equations  (4)  through  (6)  are  coupled  with  Eqs.  (1) 
through  (3)  and  integrated  simultaneously  along  a  streamline. 
The  fourth-order  Runge-fvutta  scheme  is  used  in  the  numerical 
integration.  To  account  for  the  crossflow  derivatives,  the  in¬ 
tegration  is  performed  along  at  least  two  adjacent  streamlines 
simultaneously.  The  differences  in  the  resulting  variables  are 
then  evaluated  as  the  crossflow  derivatives  and  immediately  led 
into  the  source  terms  of  Eqs.  (4)  through  (6);  see  Appendix. 
This  whole  procedure  has  been  coded  in  BASIC  language  using 
an  H-P  9836  desk-top  computer.  Calculated  flow  patterns  are 
instantly  displayed  graphically  to  facilitate  locating  the 
streamline  convergence. 

Results  and  Discussion 

Numerical  results  are  calculated  for  a  prolate  spheroid  (a/b 
=  5.91)  at  =  45  m/s  (147.6  ft/s)  and  a  -  10  and  20  deg. 
The  flow  is  largely  turbulent.  In  the  calculation,  the  flow  is 
tripped  to  turbulent  at  s/a  =  0.1,  where  the  initial  conditions 
H  =  1.35,  0(|  =  0.0005,  and  1  =  0  are  applied.  Results  of 
the  streamwise  shape  factor  H  and  the  momentum  thickness 
©n  along  streamlines  of  rj  =  1,3,  10,  and  70  for  the  case  of 
a  =  10  deg  are  shown  in  Figs.  3  and  4,  respectively.  The  ©n 
values  rise  significantly  in  all  four  streamlines  as  they  proceed 
downstream.  Nevertheless,  the  downstream  flow  exhibits 
reasonable  stability  as  indicated  by  declining  H  values  in  that 
region. 

Figures  5  and  6  show  the  calculated  results  of  the  stream- 
wise  skin  friction  coefficient  Cf,  and  the  direction  of  the  wall 
shear  stress  A,  along  the  above  streamlines.  Comparison  of  the 
present  Cft  values  with  those  obtained  by  a  finite  difference 
scheme8  and  by  experiment9  is  made.  Agreement  is  within  the 
range  of  the  data  scattering;  see  Fig.  5.  The  crossflow 
characteristics  are  represented  by  the  value  A  (A  =  tan  (3). 
Generally,  A  takes  on  a  positive  value  in  the  upstream 
region  and  turns  negative  downstream;  see  Fig.  6.  This  in¬ 
dicates  that  the  effect  of  viscosity  makes  the  streamline  steeper 
in  the  upstream  region  and  flatter  downstream  compared  with 
the  external  inviscid  streamlines  based  on  the  potential  flow 
theory;  see  Fig.  7.  The  trend  is  consistent  with  previous 
findings. I0,  1 1  The  deflection  of  streamlines  in  the  leeside  that 
causes  vortex-type  separation  is  clearly  portrayed. 

Contrary  to  intuition,  the  deviation  of  viscous  streamlines 
from  the  inviscid  pattern  does  not  stem  from  the  shear  stress 
derivative  in  Eq.  (3),  but  rather  is  attributed  to  the  deficiency 
of  the  velocity  inside  the  boundary  layer.  In  Fig.  8,  the 
crossflow  velocity  profiles  and  the  resulting  shear  stress 
derivatives  are  plotted  at  five  locations  along  the  streamline  rj 
=  3.  At  z  =  d,/d,  where  the  layer  is  set  for  the  present  com¬ 
putation,  the  numerical  values  for  3T^/3z/yM2P  from  Eq. 

(22)  are  negative  in  the  upstream  region  (Points  A,  B,  and  C  in 
Fig.  8)  then  change  to  positive  downstream  (Points  D  and  E). 
In  fact,  these  friction  values  provide  resistance  to  the  change 
of  the  streamline  direction  from  its  external  pattern. 

The  change  in  the  sign  of  the  wall  shear  angle  from 
positive  to  negative  implies  that  a  zero  crossflow  skin  friction 
line  exists,  since  Cfr|  =  ACf,.  The  calculated  crossflow  skin 
friction  coefficients  along  these  lines  are  presented  in  Fig.  9. 
Based  on  the  framework  of  the  present  theory,  this  zero 
crossflow  skin  friction  line  (portrayed  in  Fig.  10)  is  apparently 
not  a  line  of  separation.  This  is  in  agreement  with  experimental 
observation.12 


s/2a 

Fig.  3  -  Streamwise  boundary  layer  shape  factor  ot 
flow  over  a  prolate  spheroid  at  o  =  10  deg  and  V  ,  -  45  m  s 


Fig.  4  -  Streamwise  momentum  thickness  of  flow  over  a 
prolate  spheroid  at  a  =  10  deg  and  x  45  m/s 


Fig.  5  -  Streamwise  skin  friction  coefficients  of  flow  over 
a  prolate  spheroid  at  a  =  10  deg  and  V„  =  45  m/s 


Fig.  6  -  Direction  of  wall  shear  stress  of  flow  over 
a  prolate  spheroid  at  a  =  10  deg  and  V„  =  45  m/s 
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-  CALCULATED  STREAMLINES 

-  -  EXTERNAL  INVISCID  STREAMLINES 
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Fig.  7  -  Comparison  of  calculated  and  external  imiscid 
streamlines 


1  Crossflow  velocity  profile,  v/U 

2  Crossflow  shear  stress  derivative. 
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Fig.  8  -  Crossflow  velocity  profiles  and  resulting  shear  stress  derivatives  along  streamline 
rj  =  3  of  flow  over  a  prolate  spheroid  at  a  =  10  deg  and  =  45  nv  s 
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Fig.  9  -  Crossflow  skin  friction  coefficients  of  flow  over  a 
prolate  spheroid  at  a  =  10  deg  and  V„,  -  45  m  s 
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The  calculated  moral  I  flow  pattern  it  shown  in  I  ip.  10. 

The  flow  separation  it  detected  near  the  rear  stagnation  region. 
The  retult  it  compared  with  that  obtained  by  the  limiting 
streamline  pattern  bated  on  the  measured  wall  thear  tirestet.1 ' 
Agreement  it  very  good.  At  the  line  of  teparation,  which 
retultt  front  the  formation  of  the  envelope  of  viscous 
streamlines,  the  streamlines  from  the  windside  yield  highei 
stream  wise  wall  friction  values  than  those  from  the  leeside. 

This  implies  that  the  windside  streamlines  possess  more  energy 
and.  therefore,  will  roll  over  those  from  the  leeside.  A  similar 
conclusion  was  reached  from  the  analysis  based  on  the  dif¬ 
ference  in  velocity  components  at  the  edge  of  the  boundary 
layer  normal  to  the  separation  line." 

The  result  of  the  How  pattern  at  a  =  20  deg  is  shown  in 
Fig.  1 1.  The  line  of  separation  emerges  further  in  the  upstream 
creating  a  larger  separated  region  in  the  leeside  of  the  How. 
Although  a  direct  comparison  between  calculated  and 
experimental  results  cannot  be  made  because  of  a  lack  of 
published  data  for  this  particular  angle  of  attack,  the  present 
result  does  fall  in  between  the  measurements  of  the  a  --  10-  to 
30-deg  cases. i: 
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Fig.  10  -  Flow  separation  pattern  on  a  prolate  spheroid 
at  a  -  10  deg  and  V„  =  45  m/s 
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Fig.  II  -  Flow  separation  pattern  on  a  prolate  spheroid 
at  o  -  20  deg  and  =  45  m/s 
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A  three-dimensional  integtal  boundary -layer  approach  i- 
developed  and  coupled  vvnli  the  streamline  method  tot 
theoretically  determining  the  vortex-type  How  sepaiation.  1  he 
deflection  of  leeside  streamlines  that  causes  the  vortex  type 
separation  is  found  to  be  attributed  to  the  deficiency  ot  veloci¬ 
ty  inside  the  boundary  layer.  The  calculated  theoretical  results 
compare  very  well  with  the  experimental  data  foi  the  case  ot 
the  prolate  spheroid  at  a  10  deg.  For  highei  angle  of  attack, 
the  approach  yields  qualitative  agreement.  The  /eio  crossflow 
skin  friction  line  exists,  but  it  is  not  a  line  of  separation. 
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Appendix 

Functions  tor  Cl,.  and  ij, 

The  detailed  expressions  tor  functions  (.  .  and  to, 
Eqs.  (13)  through  (15)  are  listed  below. 
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Errata:  An  Integral  Prediction  Method  for 
Three-Dimensional  Flow  Separation 
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Equation  (12)  on  page  3  of  the  paper  should  read: 

,  3h„ 


K2  = 
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yM2P  \ 3x 


cos  Q  +  ~  sin  0  j 
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+  f 

dx 

Equations  for  B]^  and  Q22  in  the  Appendix  of  the  paper  should  read 
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